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The estimates of Rayleigh fading channels are rapidly changing in multi 
carrier direct sequence code division multiple access (MC-DS-CDMA) 
multiplexing systems. The most widely accepted answer to this issue is the 
conventional solution least square (LS) or mean square error estimator 
(MMSE) using the recursive least squares algorithm (RLS) or the least mean 
squares (LMS) algorithm. In much of the previous work, only one Kalman 
filter was used for estimation. In this paper, a Kalman filter is used with a 
Madgwicks filter together to satisfy the fading problems. However, this 
requires a priori evaluation of auto regressive (AR) parameters. A standard 
solution involves the first matching of the auto-completion function of the 
applying the AR method to Jakes' problem and then tackling it (YWE). Even 
more the results procedure is limited to crowd constraints and is related to an 
AR+ process of noise, an approximation considered. In fact, depending on 
simulation findings, high-AR models outperform conventional models on 
the basis of spectral estimate and bit error margins (BER). Nevertheless, in 


order to save costs of computing, the 5-D model of AR is a possibility. The 
proposed process outperforms edge of art competitors in terms of bit error 
rate as demonstrated by results. 
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1. INTRODUCTION 

Multiple access by multiplying code division is an effective synthesis of multiplicity with frequency 
division in orthogonal directions and widespread data symbols in the frequency domain [1]-[3]. It achieves 
high data rates and high overcurrent capability with extracting quasi-optimal multiplicity of multi-paths, 
avoiding inter-symbol interference, and efficient bandwidth utilization, also, the use of high-permeability 
filtering techniques, and the implementation of spatial multi-tails [4] will upgrade the communication link via 
space-frequency, space-time or space-time-frequency. Therefore, it can be predicted that the combination of 
these two techniques is one of the main options in wireless telecommunications [5]. 

The receiver design usually requires information about the channel mode [6]. According to the multi 
carrier direct sequence code division multiple access (MC-DS-CDMA) ordinary receivers [7], whose 
includes a congruence along with each carrier and subsequently a combination of the largest ratio maximum 
ratio combining (MRC) [8], [9]. In addition, blind techniques need more viewing windows and greater 
complexity than instructional techniques so, variable-time fade as a result of the Jakes model, channels are 
often depicted as random processes with a restricted doppler power spectrum [10]. In directly estimating 
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fading processes they use a least recursive squared solution or the least square algorithm (LMS) and least 
squares recursive algorithm (RLS) [11], [12]. 

Kalman filter is an optimised recursive filter used in many areas for estimate optimal state of core 
variable [13], [14]. Commonly, to explain the development of a signal, Kalman filtering using an auto- 
regressive (AR) modelling is used of faded process time shows that BER is superior to the LMS independent 
of the model and RLS based on. The coefficients and order rank of the model must be estimated to model the 
data. Additionally, the AR ordered modeling selection should be examined. There are four common methods 
for estimating the coefficients of this parameter model: Yule-Walker (YW), Burg, least mean squares, and 
least squares forward (LSF) [15], [16]. 

On another hand, the strictures of AR can also be assessed by receiver's signal. Between the 
available techniques, Tsatsanis et al. [17] have suggested the estimation of AR Characteristics from the 
performance estimation perform of the channel data covariance, using a standard yol-walker estimator. 
Nevertheless, the estimation of AR characteristics fading. Gao et al. [18], the parameters of the AR model are 
estimated using the maximizing expectancy (EM) algorithm using the Kalman harmony. Jamoos et al. [19], 
two-element Kalman filter-based structures are proposed to estimate the combined MC-DS-CDMA channels 
that are vanishing and their respective AR characteristics [20]. 

The fading process is estimated alongside any carrier using a Kalman filter, while the AR 
specifications of the vanishing process are estimated using a second estimate. Indexing and abstracting 
services rely on the title's correctness, extracting keywords for cross-referencing and computer search [21]. 
A document with an inappropriate title will very certainly never reach the target audience, so be precise. 


2. MC-DS-CDMA SCHEME 
If N is the data, the AR process introduced in x(1)...x(N) is available then [3]: 


Xn = [X(N — 1),...,x(N - n)] (1) 
En = [E(N — 1),...,E(N —n)] (2) 
E(N —k) = [e(N — k), ..,e(N — k —m+1)]" (3) 
x(N — k) = [x(N — Kk), ...,x(N — k — m + 1)]" (4) 


where En, Xn are matrices of m Xx m. For k = 1,2,...,n if (n =p) and (m2 2n), then n — p the first 
column of the matrix B = [X„ |En] depends linearly on E,, and to p the last column of X,,. Thus, B has n — p 
a single value of zero and a non-zero value. Suppose p1 = fz 2 + = B, and a, = az È = = a, represent 
the values of B and X, respectively. If (using the perturbation theory) one can obtain 6 values of a single 
value, one can obtain a threshold value a, = 6 > a1, then the effective rank X,, is equal to k. Therefore, B 
is: 


B = A +E = [X,|0] + [0|E,] (5) 
using disturbances theory and relation (5), the following relation exists between E, and A: 

|Bp+1 Z pail < & = IE (6) 

Ap+1 S &1 + Boss (7) 


if æp > 6 = & + Byy1, then a,,, < 6 < a, will be. In this case, the effective order p is obtained by the 
matrix X,,. Such a solution for 6 requires knowing £; and f,,,,. Using the theory of statistical tests and using 
soft domains, 2 matrices for a level with semanticsa, €, are limited between €,, €,. It is shown in [20] that 
with the probability 1 — a the lower limit €, and the upper limit of ¢€,, are given as [21], [22]: 


T < £1 S Vnc(m)o, = Eu (8) 


0<1t<VJc(m)o, = E, (9) 
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the c(m) represents the value of the semantic level c(m) of a distribution x? with m degrees of freedom. The 
value of 6,41 is unknown in general. In order to avoid the problem of estimating 6,41, it is assumed that p + 
1 = n and B’ is defined as: 


B' =[x(N-1)-e(N-1)  x(N—2),...,x(N-p-JD] (10) 


the first column of B” is a linear combination of other columns, so B’ has p and 6,4, = 0. On the other hand, 
we can write B’ as B' = A’ + E' = Xp+1 + [-E(N — 1)|0], so £ = ||E(N — 1)|l2 and will be according to 
the relationships (7) and (8): 


a'p+1 < & = [EWN - Dll < Vemo = £, (11) 


therefore, if a’, > 6 = €,, thena’, > ô 2 d&'p+ 
the relations (8) and (9), it can be seen that ¢, and €; depend on o,. To calculate the variance of the 
remaining error, the following error is first given to the a = [aæ1, @2, 0, Ay)’. The amount of variance 
remaining is calculated as in (13). 


, and the effective order of the matrix X, is obtained. Using 


X,a(n) = x(N) (12) 


mse, = ||x(N) — X,a(n)|I? > o? = = 


(13) 
However, according to Brorsen, the relationship between the input noise variance model and the 
remaining variance of a relationship is as: 
= 2 mr Ea 2 — G 
RV (q) = Oe MaA vi) > Oe IRa-vi) (14) 
that is, v; is the experimental values depend on the model estimation method. Here, the LSF method is used 


to estimate the parameters and v;is corrected in [22] for the LSF method and the new formula v; = 
waan given for it. Applying these changes causes the value of []{7(1 — vi) to be greater than the 
previous method due to the presence of the coefficients gẹ? in the denominator of the variance remaining. 


3. ASSESSMENT OF CHANNELS USING KALMAN FILTERING 
3.1. Rayleigh fading channels AR modeling 

Stochastic features of the mth procedure of carrier loss (h m (n)) rely on highest frequency of 
doppler shift: 


fa = vfc/c (15) 
when v denotes the mobile velocity, f c denotes c represents the speed, while f is the fundamental carrier 


frequency of light. PSD related with in-phase or quadrature fading, as described by [5], is band-limited and 
U-shaped theoretically. Additionally, it displays fd has two distinct highs and lows. In the following manner: 


uk 

a < 

Pant) -| nfav1-CF/fa) IFI < fa (16) 
0 else where 


ACF (normalised discrete-time autocorrelation function) therefore meets the following conditions: 


Rin (N) = Jo(2tfaTy|n1) (17) 
where fdTb is the doppler rate and J O is the first-order zero-order Bessel equation. The fading mechanism 


along the mth carrier may be represented using the AR(p) process, which has p spikes in the frequency 
domain: 
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hm (Nn) = — Lier aihn (Nn — Í) + Um (n) (18) 


where a i _(i=1,...,P) signifies the parameters of the the sophisticated white Gaussian driving process with 
homogeneity of variance _u2 all along carriers of the AR model. 


The following equation describes the link between the AR variables and the fade process ACF R 
hh(n): 


Ron (n) = — Vier GRinn-i) n21 (19) 


in (12) a matrix may be used to represent the data n = 1, ..., P as: 


Ran (0) Ran(-1) Rup + 1)] pa, Ran C1) 
Ran) Ran(0) Ran(—P + 2)))@2] _ _ Ran (2) (20) 

Ran - 1) Ran(P—2) > Ra) Jel Raa) 
og = Ryn (0) + Dir Gi Ran(-D (21) 


The AR parameters may be determined from R hh (n) in (21) by resolving the YWE, 
i.e. a=-R hh(-1) v, where R hh(-1) indicates the reverse fading processes autocorrelation matrices. Indeed, the 
YWE is ill-conditioned for all but the smallest AR model orders owing to the band-limited character of the 
doppler fading processes spectrum. As a result, and for simplicity's sake, prior investigations (e.g., [8], [20]- 
[22]) concentrated on only AR designs with low order. As indicated in part I, Baddour et al. [23] proposed a 
straightforward heuristic for reducing the underlying criterion of R hh by introducing a very modest a 
favourable attitude toward its major diagonal. Thus, the resultant AR (p) process's initial p+I autocorrelation 
lags meet as (22). 


Ryn (0) +e n = 0 


Ran (n) n = 1,2, ...,p (22) 


Rn (n) = { 


The magnitude of the additional bias is mostly determined by the doppler ratio f d T b, as illustrated 
in [23]. Table 1 contains representative values of variable and a trade-off among both enhanced restriction 
number of R (hh) and the model's bias. 


Table 1. Typical bias doppler rates with extra bias 


Variable Speed (rpm) 
0.001 10° 
0.005 10° 
0.01 107 
0.05 108 

0.1 10° 
0.2 10° 


3.2. Detection and prediction of discoloration 
As in (20), the objective is to calculate the fading rate of the mth carrier sequences h m (n) 
approximated through a pth order AR process. p1 state vector is specified as (23): 


h(n) = [h(n) hín- 1) .. h-pt+1)]7 (23) 


it is shown the carriers subscript is omitted for the purpose of clarification and simplicity of presenting, and 
the formula (23) is provided by: 


h(n + 1) = Oh(n) + gu(n) (24) 
where the following definitions apply to the transition matrix, the vector g, and the noticing vector h(n) are 


as: 
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Kalman filtering may be used to estimate the state vector h(n+1/n) given the collection of data provided 
below and Its variance is then defined by (27): 


a(n) = y(n) -b MÂ- 1) (26) 
C(n) = Ela(n)a*(n)] = bP C = 1) b(n) +0? (27) 


where P(n/n — 1) denotes to a at time n, the priori error covariance matrix. According to this formula, the 
Kalman profit is computed: 


C(n) = oP [E — 1] bec) (28) 

K(n) = oP |- 1] bcn) (29) 
the estimation of the state vector h (n + 1/n) and the fading process h (n + 1/n) are: 

h (n++) = 0h (n+=)+K@an) (30) 

ñ (n+~) = 0h (n+-) (31) 


h(O/-1)=0 and P(O/-1)=IP are the beginning values in state vector and error covariance template, 
respectively, in this case. Madgwick filter is an estimator based on the inputs. So that, multi-axis signal is 
estimated based on quaternion representation as follows [24], [25]: 


Ut = Ede-1 + GGrAt (32) 


VF (EanXs,) 


by = ae ead SES 
Pe yp EatXs, dil 


1 
Git = 564-1 @ Xu, (33) 
in (12), Eqs and &q;_, are the orientation in global coordinates, G, relative to the network’s local coordinates, 
L, at time t and t — 1. Also, &q, is the rate of change in signals. In (33), the signal Xo, = [0 Wy Wy @,] 
is the rotation angle that is given by the signal and the signal Xs, = [0 a, a, 4,] is the signal along axis 


x, y, and Z, that is given by the accelerometer with quaternion description at time t. The parameter p; is 
shown in (34). 


= La 
He = (@ + Be) [lédee| At (34) 
Where qœ, and a; are the zero-mean white Gaussian noise in the signal. In cases where the signal is 
estimate, alternative vectors X,, = X — €, and X, = X, — e, can be used instead. Substituting the error e,, 
and e, as measurement noise into Sẹ and sẹ gives new models as in (35) and (36). Here, the threshold of 
reconstruction error 6,, and ôg are used as as in (37). 


PREIE 7 o _. _VF(Gan%s) 

cat = 5 GIt-1 B Xu, — Me ranks II (35) 
He = (Cy + ese) [láda l| At (36) 
He = (bw + 8st) |léġol| At (37) 
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4. EXPERIMENTAL RESULTS 

In terms of BER, the MC-DS-CDMA performs in the framework is shown utilising a Kalman 
network estimate under the realistic Jakes model with high-order AR conditions with varying situations in 
which the pace of fading. Additionally, the Kalman estimator's performance is compared to that of the more 
conventional LMS and RLS estimators. The suggested solution considers a system with K=10 active users 
who each have a golden value N=3,1. The number of carriers is set to a reasonable figure (i.e. M=1 or M=3) 
in accordance with the 3rd level CDMA2000 standards [2]. Furthermore, fading steps h m (n) (m=1,...,M) are 
created using a modified Jakes model and normalised to have a unit variation, i.e. _h2=1. 

The experimental plot illustrates the expected amplitude and phase of the fade processes in front of 
very first carrier utilising AR (5) structure, the LMS predictor, and RLS interpolation method. As shown in 
the figures, the Kalman filter-based estimate outperforms the LMS and RLS estimators significantly. 
Figure 1, shows depicts the system act of the MC-DS-CDMA as a function of SNR, when Kalman method 
path estimator achieves significantlyless BER than the LMS and RLS-based estimators, which tend to a fault 
floor. 


—e— Proposed Methdos o.oo 
------ Cross Kalman 


BER*10%-3 
> 


=- ©. ~ Madgwicks Estimate 
—e— Kalman Estimate 
—— Known Channel 


Figure 1. For a variety of carriers, the result of the MC-DS-CDMaA structure of multiple channel estimators 
vs the doppler rate is shown. M is equal to one and M is equal to three. SNR=5 decibels 


By leveraging AR models to take use of the channel statistics, our proposed channel estimator 
outperforms the design LMS and RLS-based estimators. Increase the amount of carrier from M=1 to M=3, 
and you get a large increase in frequency diversity as a bonus. According to the Figures 2-5, the BER is 
affected by the different channel estimators when there is varying fading rates. As f d T b0.05, the Kalman 
estimator outperforms the others, particularly for high-order AR model. Figures 4 and 5 show that a model 
with AR (20) has lower BER than one with AR (5), making it the better choice. 


—¢— Proposed Methdos 


- -@ - Cross Kalman 

~ ©. ~ Madgwicks Estimate 
—e— Kalman Estimate 
—l— Known Channel 


BER*10%-3 


Figure 2. The performance of the proposed algorithm comparison of the MC-DS-CDMA layout with various 
channel estimation methods to the transmitor doppler level M=1 and M=3 is shown against the doppler rate. 
SNR=15 decibels 
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—¢e— Proposed Methdos 


- — - Cross Kalman 

~ ©. ~ Madgwicks Estimate 
—e— Kalman Estimate 
—m— Known Channel 


Figure 3. BER efficiency of MC-DS-CDMA scheme using different channel estimators as a parameter of 
doppler ratio for M=1 and M=3 carriers. SNR=25 decibels 


Proposed Methdos 
- = - Cross Kalman 
--&. — Madgwicks Estimate 
—e— Kalman Estimate 

50 60 
—l— Known Channel Dopper Rate 


70 100 


Figure 4. BER efficiency of MC-DS-CDMA scheme using different channel estimators when the quantity of 


transporter is M=1 or M=3 


—¢— Proposed Methdos 


- -® - Cross Kalman as 

~ ©. ~ Madgwicks Estimate Sc O 

—e— Kalman Estimate pice ee’ 
—m— Known Channel 70 80 90 100 


50 60 
Dopper Rate 


Figure 5. Channels M=1 and M=3, BER efficiency of a MC-DS-CDMA regulation with different duct 


estimators f d T b=0.01 


Multiuser QS-CDMA signal detection technique using the QRD-M algorithm and the kalman 
algorithm to monitor the latency and channel coefficients of the individual users. In general, the BER and 
ASE of kalman channel/delay estimation methods driven by MMSE and decorrelator-based judgments are 
lower when QRD-M multiuser detectors are used with the kalman method. The QRD-M method beat the 
MMSE and decorrelator detectors in simulations, even with a minimal number of pathways (relative to a 


Indonesian J Elec Eng & Comp Sci, Vol. 27, No. 1, July 2022: 262-270 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 o 269 


complete search). However, the BER does not alter the algorithm's outcomes, but only how many times the 
kalman filter is alerted that a value is in the hash table without it. 


5. CONCLUSION 

The MC-DS-CDMA systems and etiolation channels, this work investigates the speculation and 
equalisation of pilot symbols. Using two Kalman filters, one on each subcarrier, and a Madgwicks filter 
frame, to predect fading process and associated AR factors at pilot symbol points. A variety of methods are 
employed in this study to assess the suggested structure's performance. The hypothesized dual Kalman filter 
and Madgwicks content estimate beats the LMS and RLS shows in simulated results. It has also been shown 
that low-pass interpolation beats spline and linear extrapolation. 
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